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ABSTRACT 

Free vibration of laminated composite beams is studied. The effect of interply delaminations on 
natural frequencies and mode shapes is evaluated both analytically and experimentally. A generalized 
variational principle is used to formulate the equation of motion and associated boundary conditions for 
the free vibration of a composite beam with a delamination of arbitrary size and location. The effect of 
coupling between longitudinal vibration and bending vibration is considered. This coupling effect is 
shown to significantly affect the calculated natural frequencies and mode shapes of the delaminated beam. 

NOMENCLATURE 

A beam cross-section area 

b half breadth of rectangular beam 

d half depth of rectangular beam 

d l delamination location along the thickness of the beam 

E composite bending modulus 

T 

f nodal displacement vector, f = (w,^,u) 

G composite shear modulus 

h laminate thickness 
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I cross-sectional area moment of inertia 

J functional integral 

[K ] global stiffness matrix 

[kj] element stiffness matrix 

{ length of beam 

left crack tip location 
C 2 right crack tip location 

tg length of subbeam (= f/M) 

M number of subbeams 

[M ] global mass matrix 

[m>] element mass matrix 

t ply thickness 

u- displacement components (= u, v, w) 

u(x) longitudinal deflection amplitude 

u nodal longitudinal displacement vector 

V volume 

w(x,t) bending deflection 

w(x) bending deflection amplitude 

w nodal transverse displacement vector 

/? shear angle 

$ shear angle amplitude 

ji nodal shear angle vector 

6- Kronecker’s delta, = 1 for i = j and =0 for i * j 

p mass density 

e- strain tensor component 


2 


ay stress tensor component 

(aj c natural frequency of delaminated beam 

u; uc natural frequency of undelaminated beam 

£ =x/t 

7 position coordinate along a subbeam 

= d/ at 

' = 3/dx 

INTRODUCTION 

The use of composite materials has been increasing consistently in aerospace applications due to 
their high specific stiffness and strength. However, the mechanical properties of composite materials may 
degrade severely in the presence of damage. This damage may occur as a consequence of imperfections 
introduced during the manufacturing process or it may result from external loads occurring during the 
operational life, such as impact by foreign objects. 

A considerable amount of research effort has been expended in attempting to understand and 
model the failure processes of composite materials. Such a task is complicated by the fact that damage in 
composite materials may grow as a combination of several failure modes such as matrix cracking, fiber 
pullout, fiber fracture, or delamination . 1 Delamination is one of the most common failure modes of 
composite materials, and can potentially be the most damaging. Numerous researchers 2 12 have studied 
various aspects of the delamination process analytically and experimentally. Most of these studies were 
focused on either the initiation or the propagation of delamination under in-plane static or fatigue 
loadings. In comparision, little research has been done to determine the changes in dynamic response of 
composite materials caused by delamination. 

Recently, several studies have been conducted in the area of vibration behavior. The effect of 

1 3 

delamination on vibration properties of composite beams was studied by Grady and Meyn and Lee 
et al . 14 Their experimental measurements showed that the natural frequencies and damping of composite 
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test specimens were significantly affected by the introduction of an interfacial delamination. Due to the 
difficulty of mathematical modeling, the results of these tests were compared to calculations from finite 
element approach only. The possibility of determining the integrity of composite beams from measured 
dynamic behavior was also proposed. 14 Free vibration of a laminated composite beam with a one- 
dimensional delamination was also studied by Ramkumar et al. 15 Their analytical model predicted a 
severe drop in the natural frequencies due to delamination, which was found to be inconsistent with 
experimental results. They attributed this discrepancy to the effect of contact between the delamination 
surfaces. However, this explanation has never been validated. Recently, Wang et al. 16 and Yin and 
Jane have improved this analytical model by including the effect of coupling between the longitudinal 
and flexural motion in the delaminated plies. With the inclusion of coupling, the results were in good 
agreement with experimental measurements. 

The objective of the proposed work is to further develop an effective mathematical model for 
determining the natural frequencies and corresponding mode shapes for composite beams with a 
delamination. This model is based on cracked beam theories 18 " 20 derived using Timoshenko beam theory 
and general kinematic assumptions which account for the coupling effects mentioned above. A 
corresponding experimental study is also carried out to verify the analytical predictions. The use of this 
new theoretical model to identify the delamination position and size from changes in the natural 
frequencies and mode shapes of the beams is also discussed. 

THEORETICAL DEVELOPMENT 

The geometry of the delaminated beam is shown in Fig. 1. The delamination is assumed to be 
uniform through the width of the beam, and located arbitrarily, as defined by the parameters f 1? C 2 > an d 
d|. The objective of the present theoretical model is to relate these parameters to the natural frequencies 
and mode shapes of the delaminated composite. For clarity, the theoretical formulations for the 
delaminated and undelaminated regions are derived separately, as follows. 
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Undelaminated Region 

In the undelaminated regions, indicated by beams 1 and 3 in Fig. 1, only transverse vibration is 
considered. The assumptions using the Timoshenko beam theory to include shear deformation effects are 
summarized as follows: 


u y = °> u i = w(x,t), u x = -z(w' + P) 
Px = °> Py = °> Pz = P ( x > fc ) 


e xx 

= 

u', 

ff xz - -0 

Syy 

= 

e zz 

= 

V€ 

^ C XX 

e xy 

= 


= 

0 

a xx 

= 

^xx 

( x > 

Z,t), <7 — <7 xz (x,Z,t) 

ffyy 

= 

°zz 

rr 

II 

d 

V! 

est 

H 

O 

X x 

= 

x y 

= 

X 

N 

II 

O 


where Uj are the displacements referring to Cartesian axes x, y, z; a and e represent stress and strain, 
P is the angle due to shear force, and X ; and pj are the body forces and velocity components, 
respectively. 

In this work, the Hu-Washizu principle is modified to include the virtual work done by the inertial 
forces. This yields the following functional: 



i, 


/jppi - i /0PiP| - A(ey) + 
2 




ky + XjUj 


d v + J s gjUjdS + J s gi ( Ui - u ; )dS 
*^1 2 


f dt 

( 2 ) 


where a- and e- are stress and strain components, and Xj and pj are the body forces and velocity 
components, respectively. (The indices ij — 1,2,3 refer to the orthogonal directions x,y,z, respectively.) 
p is the mass density, A(ejj) is the strain energy density function, the gj are the respective surface 
tractions, V is the total volume of the system, S is its external surface, Sj is the traction prescribed 


5 


boundary, and Sg is the displacement prescribed boundary. The overbarred quantities gj and Uj 

denote the prescribed values of surface tractions and surface displacements, respectively. 

The functional J of Eq. (2) has stationary values for the actual solution for the independent 
quantities Uj, Pj, £y, and Therefore, from variational principle, for arbitrary independent variations 
of (within conditions ^u(t^) = = 0), £p-, 6e- t1 and 6a ^ the first variation of the functional J 

vanishes, i.e., S3 ~ 0, and is listed as follows: 


sj . r * 

h, 




+ x i " PPi)* u i + (*ij - 





/ 1 


5 ij - jfoj + u j,i) 

6a i ) + 

p - 

^Pipi) pi 

^Pi 


(3) 


dV 


+ J g (gj - gj)^UjdS + J s (Uj - Uj) ^gjdS 


L dt = 0 


The kinematic assumptions (Eqs. (1)) are substituted into the above formation (Eq. (3)), whereby 
the problem is reduced to a form corresponding to the beam model. The equations of motion can be 
derived in terms of displacement w and shear angle fi by subsequently integrating by parts each term 
of Eq. (7) and then integrating over the section. This gives 

E*I(w'"' + p“') + pAw = 0 

E*I(w'" + fi") - AkG = 0 

along with the associated boundary conditions. Here, according to Ref. 21, /c(= 0.7715) is the shear 
correction factor, and E , G are the effective beam bending and shear moduli, respectively. 
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Delaminated Region 

In Fig. 1, the delaminated region is divided into the upper plies (beam 2) and the lower plies 
(beam 4). The governing equations of motion are the same as the equations (Eqs. (4) and (5)) for the 
undelaminated region. 

However, the effect of in-plane displacement u 2 and u 4 on the location of the neutral planes of 
beams 2 and 4, respectively, should be considered in order to achieve the matching conditions, 

u 2 = e 2( w 2 + 02)» u 4 = e 4( w 4 + £4) ^ 

at the left and right crack tips, respectively. The distance between neutral planes of the upper plies 
(beam 2) and undelaminated region is designated e 2 . Similarly, e 4 is the distance between the neutral 
plane of the lower plies (beam 4) and undelaminated region. In the delaminated region, longitudinal 
motion is thereby coupled with bending motion, and the longitudinal vibration in the delaminated region 
changes the bending vibration of the delaminated beam. In other words, in order to describe the bending 
motion of the delamination region (beams 2 and 4), in addition to analyze Eqs. (4) and (5), the governing 
equation for the longitudinal displacement u, 

E + Au - pAn - 0 ^ 

must also be included. 

Two mathematical models are used in the present study. In the first model (model A), the crack 
interfaces are assumed to be in contact along the delaminated region throughout the vibration, and the 
coupling effect is accounted for in both the upper and lower plies of the delaminated region. Therefore, 
coupling between in-plane and transverse displacements described occurs throughout the entire 
delaminated region. In the second model (model B), contact between the delamination surfaces is 
neglected. Therefore, in-plane motion is not affected by bending, and coupling only occurs at the 
boundaries of the delamination (i.e., the crack tips). In addition, the effects associated with impact, 
friction, and penetration between the two crack faces are neglected for simplicity. 
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APPLICATION TO CANTILEVERED COMPOSITE BEAMS 


To test the accuracy of the present analytical model, the modes of free vibration of laminated 
composite beams with uniform one-dimensional delaminations have been examined. Following the 
diagram in Fig. 1, the example considered here is that of a rectangular, cantilevered beam With uniform 
cross section, of depth 2d and breadth 2b. The delamination, of length ( 2 - f 1? is located a distance 
dj from the neutral axis. Assuming simple harmonic motion, we take w(x,t) = w(x)e* Wct , 

/?(x,t) = $(x)e* Wcfc , and u(x,t) = u(x)e* Wct which give 


and 


Ei*Ii(w. + ) - w^AjWj = 0, 

E/^Wi + 0 . ) - A/cG.*3>i = 0, 



( 8 ) 

( 9 ) 


E.*AjU. - w^Ajfjj =0, i = 2,4 U®) 

Results are thereby obtained for the natural frequencies and mode shapes in terms of delamination size 
and location. 

Test Specimen Fabrication 

Seventeen laminated plates of dimension 6 by 12 in, were fabricated from T300/934 graphite/epoxy 
prepreg, supplied by ICI/Fiberite. Each laminate was laid up in an eight-ply [0/90] 2a construction and 
cured using the supplier^ recommended cure cycle. Individual test specimens, of nominal dimension 
10 by 0.5 by 0.0435 in., were cut from the laminates with a water knife that has a tolerance of ±0.005 in, 
on all dimensions. Interply delaminations were simulated by inserting a 0.001-in. thick teflon strip 
between selected plies of each laminate prior to curing. For comparison, one laminate was fabricated 
with no delamination. 

Using a pair of calipers accurate to ±0.0005 in., thickness measurements were taken of 120 
individual test specimens cut from the 17 different laminates that were fabricated for this study. 
Seventy-five percent of the specimens had thicknesses in the range 0.042 < h < 0.045 or a variation of 
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approximately ±3.5 percent from the average thickness of 0.0432 in. From Eq. (8), w varies with h 3 ^ 2 , 
so the variation in specimen thickness results in a ±4.5 percent uncertainty in the natural frequencies. 

Vibration Testing 

The apparatus shown in Fig. 2 was used to measure the resonant frequencies for the fundamental 
vibration mode of each of the test specimens. Each specimen was clamped along half of its span, as 
shown in the figure, to simulate a 5-in. long, 0.5-in. wide cantilevered beam with a single, uniform, 
through- width delamination centered at midspan. Within each laminate, a delamination of length 1, 2, 3, 
or 4 in. was embedded along one of the four-ply interfaces shown in Fig. 3. Therefore, 16 different 
specimen types were tested, each with a different delamination size/location combination. Three 
replicates of each specimen configuration were tested, to determine the measurement variability. 

An electrical resistance strain gage (type ED-DY-062AP-350) was mounted on each test specimen, 
0.5 in. from the cantilevered support, and oriented longitudinally, to measure flexural strain. A single- 
channel digital signal analyzer was used to record the transient strain response due to an impulsive force 
applied at the free end of the specimen. The first 8 sec of the free vibration response was digitally 
sampled at a rate of 5120 samples per second. The frequency response, calculated from a Fast Fourier 
Transform of the time domain data, therefore had a resolution of 0.125 Hz, which corresponds to an 
inherent measurement error of approximately 0.375 percent for the lowest mode. The fundamental 
resonant frequency was identified by locating the first peak in the frequency response curve. 

The first mode shape of each specimen type was measured photographically. Uninstrumented test 
specimens were centrally clamped to the actuator of a 100-Ib capacity shaker, which divided the 10-in. 
long specimens into two identical cantilevered beams of dimension 5 by 0.5 by 0.0435 in. The frequency 
of the sinusoidal forcing function was manually increased until the test specimen vibrated in its 
fundamental mode, as shown in Fig. 4. Because of the symmetrical displacement shape of the test 
specimen, an effectively “clamped” boundary condition is maintained at the central support point. 
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The mode shape was photographed with a black-and-white polaroid camera by choosing a shutter 
speed such that the film was exposed for approximately one cycle of motion; that is, the shutter was open 
for a time At given by 

At - 1 (11) 

w 

where w is the resonant frequency. While the shutter was open, the test specimen was illuminated by a 
strobe light that cycled on and off at a frequency fi , chosen such that 

n = Nw ( 12 ) 

where N is the approximate number of multiple exposures of the mode shape that are required. In this 
case, N = 6 was used, which gave typical results like that shown in the photograph in Fig. 5 for a test 
specimen with no delamination. Due to symmetry, only one-half of the vibrating test specimen is 
photographed. 

Local Rayleigh-Ritz Method 

Numerical integration of the mode shapes for the undelaminated beam are normally used to 
generate the mass and stiffness coefficients used in the Galerkin procedure. In this case, however, 
numerical integration is computationally intractable, because these mode shapes include hyperbolic 
functions. This makes the Galerkin procedure impractical, particularly when many modes are required. 

To circumvent this problem, a local Rayleigh-Ritz Method is used, to calculate a piecewise 
continuous fit to the deflection shape. The displacements, w(x), $(x), and u(x) are approximated by 
using cubic and linear polynomials defined over specific segments of the structure, here it is called 
subbeam. The coefficients of the polynomials are determined uniquely in terms of the displacements at 
the end points. The displacements at a point within the ith subbeam are approximated as 


w M 

= E T (r/)w., 

o 

IA 

<3 

IA 

(13) 

m 

= Wm, 

0 < Tj < (g 

(14) 
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and 


(15) 


uj( v ) = H T (*)u,, 0 < t, < t, 

where F = [F,, F 2 , F 3 , F 4 ] T and H = [H p H 2 ] T are vectors of prescribed (shape) functions of position 
and w*, and Uj are vectors of end transverse deflection and its slope, shear angle, and longitudinal 
displacement for the ith subbeam. The shape functions (Fj)j_ 1 ^ and (Hj)j_^ 2 are listed in the 

Appendix. This piecewise polynomial interpolation amounts to a finite element solution of the 
differential equations (Eqs. (8) to (10)) for the delaminated beam. In this analysis, a local Rayleigh-Ritz 
model with six shape functions, M identical subbeams, M d subbeams in the delamination region, and a 
total of 3(M - 1) + 2(M d - 1) degrees of freedom is used. 

The free vibration eigenvalue problem is expressed as 

[KJf-w’lMJf -0 < 16 


where f = (w, jS, u) T is the vector of nodal displacements, and [KJ, [MJ are the global stiffness and 
mass matrices for the entire beam. The assembly process used to obtain these matrices is symbolically 
described by 


G.[KJ,[MJ) 


M 


E 




(17) 


where fj, [kjj and [m-j are the nodal displacement, stiffness and mass matrices, respectively, for the ith 
subbeam, and the summation extends over all M + subbeams. The mass and stiffness matrices of 
the ith subbeam, in the local coordinate system are: 


KJ = Mi 


m 6 0 0 


0 nrij 0 


0 0 0 


(18) 
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and 


N = ( E i +I i) 


k 6 0 0 


0 kj k 2 

0 k 3 k 4 


(19) 


where the individual elements of the stiffness and mass matrices are given in the Appendix. 

The stiffness and mass matrices have dimension 6 by 6 for the subbeams located in the 
undelaminated region and 8 by 8 in the delaminated region. 

The eigenvalue problem, Eq. (16), is then solved for an increasing number of subbeams, M, until a 
frequency convergence test was satisfied. 


max 




< € 


( 20 ) 


where Awj^ is the change in the first frequency from the M-subbeam to the (M + l)-subbeam 
calculation, w^ is the M-subbeam estimate of the first frequency of the delaminated beam, and e is a 
small real number. For all cases presented in this paper, 100 subbeams were necessary to achieve 

convergence with e = 2.0x10 5 . 


RESULTS AND DISCUSSION 


The validity of this theoretical model is now demonstrated by applying it to a composite beam 
structure with embedded internal delaminations of varying size and at several different locations. 
Specifically, the material is a graphite fiber/epoxy matrix composite with ply properties as given in 
Table 1. The test specimens are cut from eight-ply laminates, laid up in a cross ply, [0/90] 2a 
configuration. It is important to note that, the thickness of the beam used in the following theoretical 
analysis was 0.04 in., which was determined by adding the nominal thickness (0.005 in.) of the individual 
plies. 
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Natural Frequencies 

Calculated results presented in Figs. 6 to 9 and Tables 2 to 5 show the effect of delamination size 
and interface location on the fundamental natural frequencies of the beam. Model A (solid line) is based 
on the assumption that longitudinal/flexural coupling exists along the entire length of the delamination, 
and model B (dashed line) assumes that coupling exists only at the two delamination crack tips. Data 
points indicate the average measured frequency from vibration tests of three different experimental 
specimens, which were tested as described earlier. The error bars indicate the measurement variation for 
the three tests. Data scatter increases slightly for the longer delaminations. 

The results shown in Figs. 6 to 9 and Tables 2 to 5 indicate that model A gives relatively good 
agreement between the calculated frequencies and test results at three of the four interfaces. When the 
delamination lies along the outermost interface, model B gives the best agreement with test data. This 
indicates that the longitudinal/flexural coupling in the delaminated region is insignificant when the 
delamination lies between the outermost plies of the laminate. In all other cases, coupling significantly 
increases the vibration frequency. For short delaminations, the frequencies are relatively insensitive to 
delamination size; particularly for delaminations shorter than 1 in. (20 percent of the span). Both 
analysis and test results indicate that the sensitivity (i.e., slope) of the frequencies to delamination size 
increases with the length of the delamination, for all interfacial locations. Figures 6 and 8 show excellent 
agreement between the calculated frequencies and the experimental results for delaminations along the 
midplane and interface 3, respectively. Figures 7 and 9 show that theoretical model slightly 
overestimates the measured frequencies when the delamination lies along interfaces 2 and 4. 

The theoretical results from Figs. 6 to 9 are plotted together in Fig. 10. This figure shows clearly 
that the natural frequency is most sensitive to delaminations near the midplane of the laminate, and that 
the sensitivity generally decreases as the distance between the delamination and the midplane increases. 
This could be a result of the higher degree of coupling between the longitudinal and bending motion when 
the delamination lies closer to the midplane. For delaminations very close to the outer surface of the 
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laminate (interface 4), however, the frequencies are more sensitive to delamination size, due to the lack of 
a longitudinal/flexural coupling effect, as discussed earlier. 

Mode Shapes 

The mode shapes of composite beams are also affected by the size and location of embedded 
delaminations. In this section, the calculated mode shapes are compared with those measured 
photographically for each of the 16 different combinations of delamination size/interface location 
examined in this study. 

Delaminations less than 3 in. long did not measurably affect the mode shapes. Figures 11 to 14 
show the effect of the 3-in. delaminations at each of the four-ply interfaces. If the delamination is off the 
laminate midplane (interface 1), the deformation of the laminate depends on whether the delaminated 
plies are in longitudinal tension or compression. This is illustrated by the mode shapes shown in Figs. 12 
to 14, which are unsymmetric about the line y = 0 because the delamination opens under compression, 
but not under tension. The vibration amplitude is therefore higher during the “compressive” part of the 
cyclic motion because the delaminated plies do not contribute much bending stiffness to the laminate 
when they are in compression. 

To compare the photographic results with theoretical calculations, the amplitude of the calculated 
mode shape was normalized such that the free-end displacement of the beam was the same as that 
measured experimentally. The maximum-amplitude deformed shape of each test specimen was then 
manually traced from multiple-exposure photographs of the fundamental mode shape similar to that 
shown in Fig. 5 for an undelaminated test specimen. The resulting comparisons are shown in Figs. 15 
to 18. 

The calculated mode shapes are in close agreement with the experimental measurements for the 
“compressive” part of the vibration cycle. When the delamination is along a ply interface other than the 
laminate midplane (interfaces 2 to 4), a measurable crack opening forms in the delaminated region. This 
is evident in the experimental results shown in Figs. 16 to 18, and also appears in the calculated mode 
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shapes for interface 4 (Fig. 18) although to a lesser extent. The calculations did not show a significant 
crack opening for the other interface locations. Experimental observations indicated that the relative 
magnitude of the crack opening displacement is dependent on the magnitude of the sinusoidal forcing 
function. 

CONCLUSIONS 

A formulation for the flexural motion of a composite beam containing a one-dimensional 
delamination is presented. It is based on a key kinematic assumption made to satisfy the compatibility 
requirements in the vicinity of the delaminated region. The idea is to include the coupling between 
longitudinal vibration and bending vibration in the modes of vibration. This coupling effect significantly 
increases the flexural stiffness of the laminate. 

The equation of motion and associated boundary conditions are derived under the generalized 
variational principle. The derivation procedure can be used for the cases of short beams where the effect 
of a shear stress concentration near the crack tips becomes important. The validity of the theory is 
established by examining the dynamic response of a cantilevered composite beam. The analytical 
solutions show excellent agreement with experimental results, which indicate that the flexural frequencies 
are relatively insensitive to delaminations of length less than 20 percent of the beam length. 
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APPENDIX 


The shape function for the subbeams Eqs. 18 to 20 are given by: 


F l = 1 - 3 




+ 2 




U) 
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I 2 


and 







H 2 


*1 




The elements of the mass matrix, Eq. (23), for the individual subbeams are given by: 

|m,l - JJ F T Fd,, K] - JJ H T Hd, 


and the elements of the stiffness matrix, Eq. (24), are given by: 


N ■ (E|'li)JJ B T Bd,, 

IM - E’liJJ P T Bd,, 

N - Bi’lijJ 


mi - (^fii)Jo e 7 ™. 

W - P T Pd, 

+ kQ* ^ H T Hd»f 
1 Jo — — 


where P — d/dr/H and B — d 2 /d 2 ?/F. 
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TABLE 1— PROPERTIES OF 
THE COMPOSITE BEAMS 
(EACH PLY) 



TABLE 2.— FUNDAMENTAL RESONANT FREQUENCY 
(IN HERTZ) FOR DELAMINATION 


ALONG INTERFACE 1 


«2 ' *V 

Specimen designation 

Analytical model 

in. a 

A1 

A2 

A3 

A 

B 


79.875 

79.875 

79.750 

82.042 

82.042 

'mStim 

78.376 

79.126 

77.001 

80.133 

67.363 

HH 

74.375 

75.000 

76.751 

75.285 

56.479 

3.0 

68.250 


66.375 

66.936 

47.898 

4.0 

57.623 


57.501 

57.239 

40.586 


“Delamination length. 


TABLE 3.— FUNDAMENTAL RESONANT FREQUENCY 
(IN HERTZ) FOR DELAMINATION 
ALONG INTERFACE 2 


«2 - {V 
in. a 

Specimen designation 

Analytical model 

B1 

B2 

B3 

A 

B 

M 

79.875 

79.875 

79.750 

82.042 

82.042 


79.375 

78.375 

76.626 

81.385 

68.776 

mSm 

75.126 


75.001 

78.103 

59.438 


64.001 


69.876 

71.159 

51.180 

4.0 

45.752 

49.751 

49.502 

62.121 

43.860 


“Delamination length. 
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TABLE 4.— FUNDAMENTAL RESONANT FREQUENCY 
(IN HERTZ) FOR DELAMINATION 
ALONG INTERFACE 3 


VIi’ 

in. a 

Specimen designation 

Analytical model 

Cl 

C2 

C3 

A 

B 

mm 

79.875 

79.875 

79.750 

82.042 

82.042 

■I 

79.625 

80.125 

80.625 

81.461 

75.137 

■P 

79.500 

81.875 

77.875 

79.932 

70.416 

Hfl 

75.625 

77.125 

78.125 

76.712 

65.058 

4.0 

73.376 

73.627 

70.376 

71.663 

59.131 


a Delamination length. 


TABLE 5— FUNDAMENTAL RESONANT FREQUENCY 
(IN HERTZ) FOR DELAMINATION 


ALONG INTERFACE 4 


' j|i) 
in. a 

Specimen designation 

Analytical model 

D1 

D2 

D3 

A 

B 

0.0 

79.875 

79.875 

79.750 

82.042 

82.042 

1.0 

75.375 

75.250 

77.250 

81.598 

75.834 

2.0 

69.376 

68.001 

69.375 

80.383 

71.881 

3.0 

65.375 

59.625 

N.A. 

77.698 

67.181 

4.0 

52.750 

57.876 

56.251 

73.147 

61.704 


“Delamination length. 
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z * -d 
z * d 


Figure 1 . — Geometry of a composite laminated beam with 
a one-dimensional delamination. 


Mounting 

fixture 



Figure 2. — Experimental apparatus for 
frequency measurement. 


j— Determination at ply 
\ Interface 3 



Figure 3.— Location of embedded flaw in composite 
laminate. 



Figure 4.— Experimental apparatus for mode shape 
measurement. 




(a) Photograph. 

(b) Comparison. 

Figure 5. — First bending mode shape of an undelaminated 
cantilevered beam. 
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Resonant frequency, hertz Resonant frequ ency, hertz 




Delamination length, inches 


igure 6.— First bending frequency for midplane delaml- 
tion (Interface 1). 


Figure 8. — First bending frequency for de lamination along 
interface 3. 



Delamination length, Inches 

Figure 7, — First bending frequency for delamination along 
interface 2. 



Delamination length, inches 


Figure 9. — First bending frequency for delamination along 
interface 4. 
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Interface 

1 (model a) 

2 (model a) 

3 (model a) 

4 (model b) 


_L 


0 12 3 4 

Delamlnation length, inches 

Figure 10. — First bending frequency In terms of Interface. 


Delamlnation 

length 




Figure 11 , — Measured mode shapes: effect of delamlnation 
at midplane. 


Figure 13. — Measured mode shapes: effect of delamlnation 
at interface 3. 


Delamlnation 



Normalized span position, (x/1) 


Figure 1 2. — Measured mode shapes: effect of delamlnation 
at Interface 2. 



Figure 14. — Measured mode shapes: effect of delamlnation 
at Interface 4. 
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Normalized displacement 














Normalized displacement 





(a) 1 inch deiamlnation. 

(b) 2 inch deiamlnation. 

(c) 3 inch delamination. 

(d) 4 inch delamination. 

Figure 1 7. — First bending mode shape of cantilevered 
beams with a deiamlnation at interface 3. 


(a) 1 inch delamination. 

(b) 2 inch delamination. 

(c) 3 Inch delamination. 

(d) 4 inch deiamlnation. 

Figure 1 8.— First bending mode shape of cantilevered 
beams with a deiamlnation at interface 4. 
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